Hawking radiation and thermodynamics of dynamical black 
holes in phantom dominated universe 

Khireddine Nouicer* 

Laboratory of Theoretical Physics (LPTh) and Department of Physics, Faculty of 
Sciences, University of Jijel, Bp 98 Ouled Aissa, Jijel 18000, Algeria 

Abstract 

The thermodynamic properties of dark energy-dominated universe in the presence of a 
black hole are investigated in the general case of a varying equation-of-state-parameter w{a). 
We show that all the thermodynamics quantities are regular at the phantom divide crossing, 
and particularly the temperature and the entropy of the dark fluid are always positive dehnite. 
We also study the accretion process of a phantom fluid by black holes and the conditions 
required for the validity of the generalized second law of thermodynamics. As a results we 
obtain a strictly negative chemical potential and an equation-of-state parameter w < —5/3. 
Keywords :GR black holes, dark energy theory, accretion of 

Introduction 

The discovery that the universe is currently undergoing a period of accelerating expansion, obtained 
from the observations of type la supernovae [II [2], inaugurate an exciting era of intense theoretical 
research. A variety of possible solutions to understand the mechanism driving this accelerating 
expansion have been debated during this decade including the cosmological constant, exotic matter 
and energy, modified gravity, anthropic arguments, etc. The most favored ones are dark energy 
based models and modified gravity theories such f{R) gravity [3j, and Dvali-Gabadadze-Porrati 
model of gravity [4j with an equation of state (EoS) parameter w = P/g < —1 (where p and P are 
the energy density and pressure of the cosmic fluid, respectively). The dark energy is frequently 
modeled as an homogeneous scalar field, and according to the value of the EoS parameter w, three 
different cases can be distinguished: quintessence scalar fields (—1 < w < —1/3) with a positive 
kinetic term ; cosmological constant {w = 1), where only the potential term contributes to both the 
pressure and the energy density of the field; finally, scalar fields with a negative kinetic term dubbed 
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phantom fields {w < —1). In the latter case the universe will suffer a crucial fate where the energy 
density and the scale factor diverge in finite time, ripping apart all bound systems of the universe, 
before the universe approaches the so-called Big- Rip singularity O [6]. It is also well known that 
in phantom fluid models with big rip singularity, quantum gravity effects become dominant in the 
neighboring of the big rip time [3 (H [9l [1^ . However, even phantom fluid based models suffer from 
quantum instabilities [11] and violation of the strong and dominant classical energy conditions 
|12] . the phantom fluids are favored by the cosmic microwave background experiments combined 
with large scale structure data, the Hubble parameter measurement and luminosity measurements 
of Type la supernovae |13| . 

An other important and growing field currently under investigation is related to the thermody- 
namic properties of an expanding universe jll]. Recent studies on phantom thermodynamics show 
that the entropy of the universe is negative |15| . while the generalized second law of gravitational 
thermodjTiamics (GSL) is satisfied, Sf + Sc > 0, where 5*/ is the phantom fluid entropy and 
Sc is the entropy of the cosmological horizon [T6j. An other point under debate is the influence 
of a non-zero chemical potential on the phantom thermodynamics and its relation with the GSL 
[I71[T8]. 

In this paper the thermodynamic properties of black holes immersed in dark energy-dominated 
expanding universe and the accretion process of phantom energy onto black holes are investigated. 
The first paper dealing with the later process is due to Babishev et al [19], where ignoring the 
backreaction effect of the phantom fluid on the black hole, they found that the change rate of the 
black hole mass is negative. However, in recent scenarios where the backreaction is taken into 
account 120} [21], it is found that the mass of the black hole is always an increasing function in an 
expanding Friedman-Robertson- Walker (FRW) universe. 

The organization of this paper is as follows: in Section 2 we review the exact solution recently 
obtained in [20ll2T| , and describing a black hole embedded in an expanding FRW universe. In sec- 
tion 3, we examine the Hawking radiation at the apparent horizon, and compare with magnitude 
of the phantom energy accretion process. We will show that the former is highly suppressed, par- 
ticularly at late times. In section 4, we study in an unified and general way the thermodynamics of 
cosmological black hole embedded in an expanding (FRW) universe with a general EoS parameter 
w{a), and particularly we obtain solutions realizing the crossing of the phantom divide line. In 
section 5, the stability of the solutions of section 4 under the quantum correction due to the con- 
formal anomaly is established. In section 6, we study the conditions required for the validity of the 
GSL when the black hole is immersed in phantom energy-dominated FRW universe. Particularly, 
in order to protect the GSL, we obtain a critical mass of the black hole of the order of the solar 
mass for particular values of the parameter = — /io^o/Poi where fiQ, Uq, Pq are the present day 
values of the chemical potential, the particle density and the energy density, respectively. Finally, 
we discuss and summarize our results in section 7. 



2 Cosmological expanding black hole 



The first solution of Einstein's tlieory of general relativity describing a black hole like object 
embedded in an expanding universe was introduced by McVittie in 1933 |22], and is given in 
isotropic coordinates by 



+ 2a{t)r J 

where a{t) is the scale factor and Mq is the mass of the black hole in the static case. In fact, when 
a{t) = 1, it reduces to the Schwarzschild solution. When the mass parameter is zero, the McVittie 
reduces to a spatially fiat FRW solution with the scale factor a{t). The global structure of ([1]) has 
been studied and particularly it has been shown that the solution possesses a spacelike singularity 
on the 2-sphere r = Mo/2, and cannot describe an embedded black hole in an expanding spatially 
fiat FLRW universe [23| 124] . On the other hand the McVittie solution is constrained by the non- 
accretion condition onto the central mass, and therefore is not suitable for a study of cosmic fluid 
accretion process onto black holes embedded in an expanding FRW universe. 

In the following we adopt the new solution describing a black hole embedded in a spatially fiat 
FRW universe[20l [21] 

rfs2 ^ -?^dt^ + a\t)A\r) (dr^ + r^dfi^) ^ (2) 

where A{r) = (l + ^) , B{r) = (l - a()f:) is the scale factor and Mq is the mass of the 

black hole in the static case. In fact, when a(t) = 1, the solution ([2]) reduces to the Schwarzschild 
solution, while when the mass parameter is zero, it reduces to a spatially fiat FRW solution with 
the scale factor a (if:). 

Using the areal radius r = r (l + ^^7^)^ and R = or, the metric takes the following suitable 
Painlevi^oe-GuUstrand form 




ds' = - 
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where H = d/a is the Hubble parameter and over dot stands for derivative with respect to the 
cosmic time. The term R^H^ plays the role of variable cosmological constant. Now, we introduce 
the time transformation t — y t to remove the dtdR term 



dt = F-^ {t, R) 
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where the integrating factor F {t, R) satisfy 



OrF-' = dt 
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Substituting (jl]) into ([3]) and replacing t — t , we obtain 
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The apparent horizons (AH) are solutions of h"'^daRdbR = 0, which leads to 



K = 0, 

where we introduced the Hawking-Hayward quasi-local mass 

mnit) = Moa(t). 
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A remarkable feature of this quantity is that it is coordinate- independent, and consequently is 
recognized as the physically relevant mass of the black hole. Obviously, it is always increasing in 
an expanding universe [21]. Therefore, the calculation of the change rate of the black hole mass 
will lead to opposite conclusions to that of Babishev et al. [19j. 

Discarding the unphysical branch with the lower sign in [TJ the AH are given by 



^2H " - 8GmH(t)) , Rc = ^{l + " 8Gm^,(t)) , (9) 

where Rc and Rb are the cosmological and the black hole AH, respectively. Note that the AH 
coincide at a time t^, defined by d{t^) = 1/8GMq. This coincidence takes place in a future or past 
universe depending on the kind of fluid accretion onto the black hole. 

Let us now consider that the fluid is described by a imperfect fluid with the stress-energy tensor 
given by 



(10) 



where u^' = (^,0,0,0) is the fluid four velocity and = (0,g, 0,0) is a spatial vector fleld 
describing the radial energy current. Written in terms of the comoving AH, the solutions of the 
Einstein equations of motion are |21] . 



H =^Rap, 3H+ — = -StiGRaP. (11) 

Assuming a radial heat inflow (g < 0), the change rate of the black hole mass is then 

mn = GaB^A\q\, (12) 

where A = j J dOdip^/g^ = Anr'^a'^A'^ . This relation shows that the Hawking-Hayward quasi-local 
mass is always increasing. 



3 Hawking radiation of apparent horizon 

For the apparent horizon to exist in an expanding universe we have the following condition 

This condition requires us to consider massive objects for which GtuhH is a small quantity. Using 
the present day value of the Hubble parameter, we estimate the critical BH mass as m//(t^,) < 
10^^ Mq. This is the condition for which the engulflng of the universe by the BH is prevented |25| . 
Using the areal radius R = ar (^1 + ^^7^)^, and solving for r , we obtain one physical solution given 
by 

= ^ (^Ra - GniH + ,J R\ - 2RAGmH^ . (14) 
For the solution to be real, we have to impose the condition 

Ra > 2GmH, (15) 

which is true only in an expanding universe due to the relation Ra (1 — HRa) = 2GmH- We note 
that the equality sign has been discarded in (ITSl) . since it leads to HRa = 0, which does not hold 
for rriH 7^ 0. The relation (fT5|) means that we are dealing with systems which are small compared 
to the cosmological curvature, and that distances below 2lpi are naturally excluded. 

Assuming now the EoS, p = wp, and substituting ( ITSl) in Friedmann equations, one flnds the 
following upper bounds on the densities 

p(t)< T-^, bl < (16) 

1287rm^(t)' USnmj^ {t) ^ ' 



Therefore, in case of phantom energy driven expansion, the finite increase of densities with time, 
will avoid the big rip singularity. On the other, since the BH mass increases with time, the 
upper bounds (fT6|) . become very small. As an estimate of the energy density we have p{t) < 
0.02 X 10-1° (GeV)^ , for the smallest super-massive BH detected in the dwarf Seyfert 1 galaxy POX 
52 with niH ~ 10^ M© pj. On the other hand if mn {Q ~ Mq, we have p (t) < 0.022 (GeV)^ . 
Even if we take rriH of the order of the Planck mass we have p (t) < 

We now consider the Hawking radiation from the cosmological AH. The temperature on the 
AH is defined hj Ta = \k,a\ /27r, where ka is the surface gravity. In dynamical spacetime there is 
no timelike Killing vector, and the usual definition of the surface gravity may be modified. In this 
case case the surface gravity is related to the so called trapping horizon. Here we follow the work 
of S. A. Hayward where the surface gravity is defined by K^'V[aKi,] = uKa, {a, 6 = 0, 1), with 
K"^ = —e°'''VbR the Kodama vector corresponding to the background given by Eq.Q. Evaluating 
all the quantities on the trapping horizon, we get the simplified form, ka = jDhR, where the 
metric hab is defined by ds"^ = habdx"'dx^ + R{x)dVl? . Performing the calculation, one finally finds 



(17) 



Writing itih is terms of the AH, the temperature at the AH becomes 
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Expanding f|T8|) to first order in Gitlh at the black hole AH and cosmic AH, we obtain 
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respectively. Using (|T5|) in the first factor, along with Rc < 1/if in the second factor, and 



H 



[1 + w)^ one obtains an instantaneous upper bound for temperature at cosmic AH 
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The relation (fT9|) shows that the accretion process tends to lower the Hawking temperature at 
the cosmic AH. We note that in (1201) we have not used the absolute value, as a consequence of 
the equivalence principle [28j. A similar behavior is also observed if one anticipates and use the 
scale factor f l67|) derived in section 4, for a phantom dominated era. In fact, one can see an 
unusual behavior of the temperature when approaching the time t^,. As it is shown in Fig. 1, the 



temperature at the cosmic AH increases with time, reaches a maximum, then begins to decrease, 
and stops at the end time t^,. This means that the Hawking radiation at the cosmological AH will 
decrease rapidly in favor of a hugh increasing of the accretion of phantom fluid. On the hand, the 
temperature of the black hole AH starts to fall at early stage because of the increasing of the black 
hole mass by accretion of phantom energy, and at late stage begins incraesing. However, the latter 
strange behavior is mainly due to the absolute value in the definition of temperature. 

An independent derivation of (fT6l) and (12T1) can be performed by using the following sim- 
ple arguments. Imposing positivity of the {—gtt) component of the metric, and replacing the 
quasi-local mass by a density energy we have, < —gu = 1 — 2Gmi{/R — i^^Gmn/R ~ ^ ^ 
{8TiG/3)pR^ - ^_(^£g%)pK^ ■ Hence, p < (1 - RH) / (SirGR^ /3). Using R > 2GmH, one obtains 
which is consistent with ( fT6l) . If one uses the Stefan-Boltzmann law p = aT'^, 
and the quantum mechanical relation R > 1/T [29], one finds T < (S/Ana) rUpi, which of the same 
order of magnitude as f l2T]) for mfj = irtpi. Finally, we point that the extremal case corresponding 
to Tc = never occurs, since in this case we must have Rc 
condition Rc > 1/ {2H) . 
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Figure 1: Variation of temperature associated with the black hole and cosmic AH with time for 
Mo = 0.25. 

Now, neglecting the accretion of radiation in phantom energy dominated era, and taking into 
account only the semi-classical Hawking evaporation and the phantom energy accretion term, the 
differential equation for the black hole mass reads 



dm 



H 



dt 



-4nRlaT^ + ninH, 



(22) 



where cr = A^7r^/120 is the Stefan-Boltzmann constant for massless fields with effective degree of 



freedom A^. Substituting the Hawking temperature associated with the BH apparent horizon, we 
obtain 

dniH (Tiy2 |i _ I (^1 _ ^1 _ SGmnH) (1 - w)'^ 



+ mnH. 



(23) 



dt 16n^ {1- ^yl-8GmHH) 

This is a comphcated relation for tuh, whose behavior is shown in fig. 2. Since the two terms are 
in competition, there exist a transition time, the phantom time, after which the accretion process 
dominates and the BH mass increases. Consequently, the BH does not lose but gain mass due to 
hugh accretion of dark energy. Indeed, it easy to show that the maximal rate gain mass for m^/ ^ 0, 
is {dfnH/ dt)^^^ ~ m^i. It is important to note that the accretion term becomes predominant at 
earlier times for massive BH. Next, we perform the same analysis on the variation of the mass 
inside the cosmic AH. Let us define the ratio between the radiation and the accretion term 



ATiRlaT^ 



(24) 



Substituting (fT8|) . and repeating the procedure leading to ( 12T]) . one finds 



Vc {t) < 



|1 — 3w\ / rripi 
2457607r 1 



(25) 



This result clearly shows that the Hawking radiation at the cosmic AH is insignificant, even near 
the Planck scale. 
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Figure 2: Variation of total mass associated with the black hole AH (left panel) and cosmic AH 
(right panel) with time. 



Finally, let us point an other crucial feature of the model of dynamical black hole considered in 
this paper. If the expansion of the universe is driven by phantom fluid, the black hole AH expands 



while the cosmic AH shrinks as the universe expand until they meet at Rcrit = l/(2-ff) at time t^, 
solution of (fT3|) . At times t > t^,, both the AH disappear living a proper singularity, well before 
the big-rip singularity is reached. The question if this singularity is naked or located inside the 
AH, and its connection with the violation of the Cosmic Censorship Conjecture (CCH), is still 
under debates [3T| [32] . To avoid discussing this topics, which are out of the scope of the present 
paper, and the fact that the phantom driven expansion of the universe is non singular, and that 
the radiation power can be neglected in comparison to the accretion of phantom fluid onto the BH, 
particularly at later stage, we limit the analysis in the remaining sections to the interval t < t*, 
far from the big rip singularity, and then purely classical treatment will be considered. 

4 Thermodynamics with varying w 

We now consider the thermodynamics properties of the solution described in section 2, with a 
variable EoS parameter, p{a) = w (a) p (a). The particle fluid and entropy fluid currents, and 
S'^ are given by 



where n and s are the densities of particle number and entropy, respectively. The conservations 
laws, T'!^ = 0, A^':'^ = and S'^^ = 0, computed on the background given by (j2]), give the following 
set of differential equations 



m = nu", S" = su^, 
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The solutions of the above equations are 




(29) 
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where the prime stands for derivative with respect to the scale factor, and po, hq, sq are the 
present day values of the corresponding quantities assumed to be positive definite. Here we note 
the important corrections due to the presence of the black hole. 

In the pure dark energy-dominated universe, Mq = 0, we have Rb = and Rc = l/H. 
Consequently, using the equations of motion (fTTl) . we show that 



R'c (a) _ 3 
Rc {a) ~ 2 
Substituting in Eqs.( l29ll3Tl) we obtain [33j 
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Now, assuming that p = p (T, n), p = p (T, n) and using the Gibbs law 
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combined with 



and the continuity equation, we obtain 
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Calculating {^)^ from Eq. fl29|) . we get the equation governing the evolution of temperature 
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Solving for w{a) ^ —1, we obtain 
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Using this expression, we write the energy density as a function of temperature 
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Extracting the scale factor from (fT8l) . we obtain the generahzed Stefan-Boltzmann law 
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Let us now reproduce the expressions of the temperature and energy density in the pure spatially 
flat FRW universe. Indeed, substituting (l32ll in (HOll and (ITSjl . we obtain 



P {T) = Po 



T{a) (wq + I) 
To {w{a) + 1) 



ITTf 3 tu(n)-u)0 
2 iii(o) 



exp 



w a 



ao 



daw' {a) In a 



(41) 



and 
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which are exactly the relations obtained in |33) . 

Let us now scrutinize the behavior of the temperature when w{a) crosses -1. In this case, 
Eq. (l37|) becomes 
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where w^'Ip^L is the value of w{a) at the phantom divide line (PDL) and C is constant of integration. 
The solution of (l43l) is: 
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Here, we note that unlike the standard vacuum solution with T = [33], the solution considered 
in this paper allows for a vacuum solution with a non-zero temperature. This was expected, since 
this temperature is associated with the AH of the black hole in the absence of dark field. The 
zero temperature vacuum state is recovered by setting Mq = 0. On the other it is important to 
observe that the expressions of temperature and energy density are regular everywhere including 
the phantom divide crossing. This can be easily verified by substituting the solution of the AH at 
the PDL, R{a) C'e""''lpDL("-"o) Eqs. ffT8]l39l) 
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Obviously, Eg. fl45|) shows that the temperature remains positive when w{a) 
the following ansatz for the AH 
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with /a (a.) — > 1, we rewrite the temperature as 
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which shows that the temperature is always positive definite regardless the value of w{a). When 
Mq = the temperature is positive for w{a) > 1, negative for w{a) < —1 and zero for w{a) = —1 
P3]. 

Now, we calculate the chemical potential defined by the Euler relation 
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Using the expressions of p (a) , n(a), s(a) and T(a) we obtain 
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where 
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no riQ 

is the present day chemical potential. We note that in general the sign of /xq can be arbitrary, and 
consequently the sign of /i (a). 

The entropy of the universe can be derived from (H9l) and is given by 
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Using again the relations for p{a) , T (a) , /i (a) and n (a) we obtain 
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where the present day entropy is 
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Now, defining the comoving volume by 
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we derive from Eq. (l53l) the usual entropy conservation law 
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As we did for the energy density, we write the entropy in terms of temperature as 
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Let us now consider the expressions of the chemical potential and entropy at the phantom divide 
crossing. Using (3.21) we easily show that 
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Note that using relations ( 155]) and f l5^ we easily verify that the Euler relation remains valid 
at the phantom divide line. 

Finally, let us list the corresponding relations for the energy density, temperature, entropy and 
chemical potential in the case w (a) = w = const, 
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It is interesting to note that the dependence of the thermodjTiamics parameters on the temperature 
is of the same form as in standard thermodynamics of purely expanding FRW universe. 

Before ending this section, let us reconsider explicitly the avoidance of the big rip singularity 
in a phantom dominated universe. Using Eqs. fl26|6ip we obtain 



a{t) = ao 
H{t) 



-Ho (1 + w){t- ts) 



3(l+») 



where the big rip time is 



3{l + w){t- t. 



ts — tr\ — 



3Ho(l + w) 



(67) 
(68) 

(69) 



Using the constraint f|T3|) on H{t), in phantom dominated era, one finds 

l6GmH{t) 



{ts ~t)> 



3(kl - i: 



(70) 



Then using the FRW equations and (fT5|) . we find 



P< 



it - ts) 



-1 



SnG^mn (1 + w) 



W {ts - t) 



-1 



iTiG^rriH (1 + w) 



(71) 



Substituting (170|) in these relations we reproduce the relations ( IT6|) . Hence, the universe evolves 
towards a state where a, p, p and higher derivatives are finite, without ever reaching the big rip 
singularity. This final state is reached in a finite time given by 



SGHoMna 



OIVIoUq 

3Ho{l + w) 



l+3w 
3(l+m) 



(72) 



This behavior of the universe is similar to that observed in the framework of generalized uncertainty 
principle (GUP) corrected FRW universe [35] . 



5 Stability of the solution 



In the following, we focus on the regime of cosmic dynamics where the universe undergoes a phase 
of quasi-exponential expansion, such that H / H"^ ^ 1, and examine the stability of the solutions 
obtained above at the crossing of the phantom divide line when quantum effects due to conformal 
anomaly are taken into account. In general, the conformal anomaly is given by [36j 

T = h(^F +^1^11^ +h'G + h"nR, (73) 

where F = C^^xkC'"'^'' = - 2R^^R^"' + R^^xkR''''^^ is the square of the Weyl tensor and G 
the Gauss-Bonnet invariant, G = R^ — AR^^R^'^ + R^j^^xi^R^'^^'^ , and the coefficient h and h' are 

^ ^ 120 (47r)^ ^""^ + ^""'Z' + ^^""'^ ' ^^^^ 
h' = ^£^2 K + 11^1/2 + 62ni) , (75) 

where no, ni/2 and ni are the number of scalar, Dirac fermion and vector fields, respectively. 

Since we have the condition ( fT3l) . we just give the expressions of F and G to leading order 
in GtuhH as 

F ~ O {{GmnHf) , ~ O [{GniuHf) , (76) 

G - A2h'H^ + USb'H^GniH. (77) 



Now assuming that 

and using the conservation law 

one finds 



Ta = -Pa + SPa. (78) 



Pa + 3H{pa+Pa) = 0, (79) 



Then pa and Pa are given by 

Pa = --J dta^HT,, = I^L±P±. (81) 
J 3 

In the quasi-exponential expansion we can write Pa ^ Pa ^ Ta, and then we have 

Pa^Pa^ A2b'H^ + USb'H^GniH. (82) 



On the other hand the fluid density energy can be approximated by 



SttG 



3H^ 



+ 



A-nG 



(83) 



Taking the present day values of the Hubble parameter, Hq ~ 10^^^ eV, and G ~ 10^^^ (^V)" , 
it is easy to verify that the quantum correction are very small when crossing the phantom divide 
line. Even we consider BH masses at the end of the expansion of the order of ~ 10^^ Mq, the 
correction term to the standard result is insignificant. Hence, the solutions given in section 4 are 
stable under the conformal anomaly quantum corrections. 

6 Accretion of phantom fluid and constraints on the GSL 

We now consider the problem of the validity of the GSL when the effect of the back-reaction effect 
of the phantom fluid on the black hole is taken into account. We restrict our study to the scenario 
of a black hole with a small quasi-local mass immersed in phantom fluid dominated era. In this 
case only the cosmological AH contributes, and the total entropy consists essentially of the sum of 
entropy of the cosmological AH and the entropy of the phantom fluid in thermal equilibrium with 
the cosmological AH. Indeed, using (166|) we have 



The first term is the entropy of the cosmological AH and the second term is the phantom fluid 
entropy inside a comoving volume V. Taking the derivative with respect to time and using the 
following approximation 




(84) 




2Grnh 
H 



(85) 



along with the relations rhh 
Hayward quasi-local mass 



mh'd/a, d = a i + Hj,we obtain to leading order in the Hawking- 




(86) 



where we have used the continuity equation (l26l) and ( l66l) . In order to satisfy the GSL, S > 0, the 
quantity inside the square brackets must be positive definite, yielding to the condition 



2mH > 



2nH 



(l+w) \po 



1 

l + w 



^-1(1+^) 



(l+w) \ po 



1 

l+w 



271 [1-^ 



(87) 



Now substituting 



'l + w), we rewrite the constraint on the GSL as 



GniH > —3 



(1 



w 



(l + Sw) 



A _ 1+H 
^ 2 



(5+3w)(l+w) 
l+3w 



with ^ 



2lT 



We know from Eq. (fT2|) that the Hawking-Hayward quasi- local mass is an increasing function 
of time in an expanding universe, ttlh > 0. Then, imposing the positiveness of the RHS of Eq.( 
and assuming the positiveness of the entropy, we have one solution given by 



(5 + 3«;)(l + w) 



> ^ > 0, 



(89) 



(1 + 3w) 

which is satisfied for an EoS parameter in the range w < —5/3. Using the expression of Sq given 
in (15^ and defining = —^qUq/pq, an immediate consequence of the positiveness of the entropy 
is that w > —1 — Uq. Since w < —5/3, we get a lower bound estimate to the present day value of 
the chemical potential. 



(90) 



which show that the chemical potential in phantom energy dominated-era is strictly negative. Now, 
with w = —2 and the present day parameters, we obtain the following bound on the entropy 



So < 3.5 X 10-^GeV" 



(91) 



Or in terms of and using ( !90|) . we have 



< ao < 1 + 1.26 X IO^^Tq. 



(92) 



Now, if we take Tq ^ 10^^^ GeV we have ao < 10^'^. Here we would like to point that the bounds 
on the present day values of the entropy and the parameter cire of the same order of that 
obtained by Lima et al. ^37]. However, the calculation in this paper is performed by ignoring the 
back-reaction of the phantom fluid on the black hole. In fact, they used the Schwarzschild metric, 
which cannot describe the properties of black holes embedded in an expanding FRW universe. As 
we explicitly show, taking into account the back- reaction effect, even in a low density background, 
leads to a drastic constraint on the EoS parameter and rule out the approach based on zero chemical 
potential |18] . 



Now, using the obvious relation mn = rhn / H, we obtain the quasi-local black hole mass above 
which the accretion of phantom fluid onto the black hole is permitted 



[l + w) 



GH (1 + 3m;) 
Using the constraint on the GSL fl89l) . we obtain 

rriH > rnH,crit 



A — 1+H 

^ 2 



A 



GH 



(5+3w)(l+w) 
l+3w 



A — 1+H 
2 

(5 + 3w) 



(93) 



(94) 



In terms of ao and the present day quantities, the critical quasi-local mass is then 



mMcrit = -42.27 X 10 



-20 



1 + w + ao 
5 + 3w 



M, 







Tn 



(95) 



For large values of ao ~ 10^*^ and Tq ~ 10^^^ GeV, the critical quasi-local mass is approximately 
fTiH,crit ~ 10^^ Mq. This is a huge value, allowing all black holes in the universe to accrete phantom 
fluid. We note that in section 3, we have shown that rriH^crit is the maximal allowed mass which 
prevents the engulfing of the universe by the BH. On the other hand small values of ao give critical 
BH mass of the order of the solar mass. For instance, taking aQ = 2, w = —2 and To ~ 10~^^ 
GeV, we obt ain iTLf{ Qrit ~ 4.2il<fQ. Our results are similar to that obtained in | l37j . where the EoS 
parameter is restricted to values less than —1 and the back-reaction effect ignored. We have plotted 
the present day critical quasi- local mass with respect to w for different values of the parameter ao- 
Obviously, the value of the parameter ao is very important in determining a critical quasi-local 
mass of the order of the solar mass. Finally, if we set the chemical potential to zero, the critical 
mass reduces to 

3 (l + w) pqVH (T 
mH,cra- 2vr(5 + 3w7) To V^o, 

which is negative for w < —5/3. This result is expected from the relation ( l90l) . 
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Figure 3: Variation of the critical Hawking-Hayward quasi-local mass with respect to w for different 
values of the parameter ao. 



In summary, we have investigated the thermodynamical properties of black holes immersed in 
an expanding spatially fiat FRW universe, for general EoS parameter w{a). Particularly, we 
found that the temperature of the dark fluid is always positive regardless the value of the time 
varying EoS parameter, and that the instantaneous vacuum state is characterized by non-zero 
temperature, entropy and chemical potential, respectively. An other important result is that all 
the thermodynamics parameters are regular at the phantom divide crossing. We have also analyzed 
the accretion process of phantom fluid onto black holes with small Hawking-Hayward quasi-local 
mass, and the constraints on the validity of the GSL. Particularly, assuming the positiveness of the 
phantom fluid entropy, we have shown that phantom fluid may have a negative definite chemical 
potential in order to satisfy the GSL. We have also obtained, for an EoS parameter within the 
interval w < —5/3, a critical quasi-local mass of the black hole, above which the GSL is always 
protected. The present analysis, show that taking into account the back-reaction effect of the 
phantom fluid on the black hole, even in a low density background, leads naturally to positive 
temperature and negative chemical potential, and may contributes to resolve the controversy on 
the subject [T81I37]. 
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